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Abstract 



The mass of the charm quark is analyzed in the context of QCD 
finite energy sum rules using recent BESII e + e~ annihilation 
data and a large momentum expansion of the QCD correlator 
which incorporates terms to order a 2 {m 2 / 'q 2 ) 6 '. Using various 
versions of duality, we obtain the consistent result m c (m c ) = 
(1.37 ± 0.09)GeV. Our result is quite independent of the ones 
based on the inverse moment analysis. 
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1 Introduction 



Recently the BESII collaboration has presented new data on the total e + e~ 
annihilation cross section above the charm threshold Data in this en- 
ergy region are particularly relevant for the extraction of the charm quark 
mass, one of the fundamental parameters of QCD. The charm quark can 
be determined by comparing suitable positive moments of these data with 
the corresponding moments of QCD perturbation theory. This direct quark- 
hadron duality approach was originally applied to the charm region in Ref. 
0. A reanalysis along these lines seems indicated in view of the fact that 
apart from new data, enormous progress has been made in the theoretical 
calculation of the relevant QCD correlator in the region q 2 ^> m 2 . The cor- 
relator is now known to 0(a 2 ) and 0(m 12 /q 12 ) M, so that the question of 
convergence can be meaningfully discussed. There exist also a result to 0{pP s ) 
for the quartic mass correction J| which we will not consider for reasons of 
consistency. We believe that in the case of the charm quark mass the direct 
duality approach employed by us is less prone to theoretical uncertainties as 
the more popular one based on inverse moments , , , [0 > • ^ should 
be pointed out, that we use as phenomeno logical input only the new BE- 
SII data. This is because older data in this region are plagued by unknown 
systematical errors and appear to be mutually inconsistent. An alternative 
would have been to adjust the normalizations of the various data sets so that 
they agree for large q 2 with QCD [[10] . 

We will manipulate our data rather on the basis of QCD duality. Using 
suitable linear combinations of moments the emphasis of the hadronic integral 
can be shifted at will to experimental regions where the data errors are 
small. This technique, which was originally proposed in |ITJ| , allows in many 
cases more accurate prediction, and supplies in addition beautiful consistency 
checks. This will be explained in the following. 



2 Cauchy sum rule 

QCD duality means that the theoretical and phenomenological information 
is being related by means of Cauchy sum rule 

r R 1 If 

/ -1wIL(s)p(s)ds = -— f Il QCD (s)p(s)ds (1) 

Jso 7T 2m J\s\=R 

where Imfl(s) is defined in terms of the total e + e~ annihilation cross-section 
R(e + e~ -> hadrons) = 12n ^ Qjlmn(s) (2) 

flavors 
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The experimental charm physical threshold in Eq.(|I]) is taken from the J 
resonance mass 

s = (3.097Gev) 2 (3) 
We have included in the sum rule a polynomial weight function p(s) 

p(s) = £ a n s n (4) 

which makes the sum rule a linear combination of moments of Cauchy sum 
rules. The polynomial may be chosen in a suitable way to enhance or remove 
part of the phenomenological input in the calculation. 



3 QCD integral 

The two-point function Hqcd(s) is known to 0(a 2 ) as a series expansion in 
powers of m 2 / s up the sixth power, 



6 3 (m 2 Y / -sV 

I1qcd(s) = J2z2 A v( m ,f J ') — ln 7T 
i=o 7=0 V s ) \ r J 



(5) 



where the coefficients Aij(m, /x) may contain powers of mass logarithms ln(m 2 / /i 2 ). 
The convergence of the series is not seriously affected by the mass logarithms 
provided 



In ^— < In -f- 



(6) 



QCD 



This is always the case for the scales of \i (~ 5GeV) we use. 

At tree level, the first few terms in the expansion of Uqcd{s) are given 

by 



H-qcd(s) 



16tt 2 



20 4. —s o m 2 (m 2 \ 2 / nl m 2 nl -s\ 
--In— + 8 h — 4-8 In— + 8 In— + 
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We use the strong coupling constant a 2 and running MS mass renormalized 
at the scale /i. The lengthy full expression to 0(a 2 ) and 0(m 12 /q 12 ) may be 
found in Ref.||. 

There is also a small non-perturbative contribution arising from the gluon 
condensate, which is known to 0(a s ) ||12|| . The first few terms are 



n np (s) 



—GG 

71 



~ (-X - l(-X - i-X (-1 - ln^ + ln^ + 



12 
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The Cauchy-Integral in Eq.(ffl) needs the evaluation of 



2m J\s\=r 



s k { In — | ds 



(7) 



for j = 0,1,2,3 and k = —6, —5, —4, ...These integrals can be evaluated 
analytically with the result 



J(k ? -l,j) = — i Win — ds = — x 
2m J\s\=R \ \i l ) n + 1 

1 



jln J 



_ X R j(j-l) M R /tt 



jj? n + 1 



ln J " 



o j(J - 1)0' - 2 ) lni " 3 — 7 

6 (n + 1) 2 J JU 1 - 2 



R 

7 2 



and 



J(k = -l,j) = -L/ -fln-J] rf S 



See also Ref. HT3 



4 Data integral 

The BESII data for the total e + e~ annihilation cross section is represented 
in figure 1. 

The errors given by the authors of Ref. [|IJ distinguish between systematic 
and statistical errors. In our analysis this distinction is absolutely essential as 
the statistical errors average out almost completely in the integration process 
while the systematic errors prevail. The threshold of the charm continuum 
is at twice the D Q mass which corresponds to so CO nt = 13.9 GeV 2 . For 
definiteness we discuss here only the integration up to s max = 25GeV 2 , the 
maximum value measured by BESII. If we require for duality that the data 
within their errors agree with QCD perturbation theory then it is seen that 
s max may also be chosen somewhat lower. We define 

/•25 

h = Rcontinuum(s)p(s)ds (8) 

J13.9 

From this integral over the total cross section the contribution of the light 
quarks must be subtracted to isolate the pure charm contribution. As s > 
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13.9GeV 2 it is perfectly safe to use QCD perturbation theory for contribution 
of the light flavors. We use the four loop result for the massless correlator 



I 2 = 2/ (l + — + ..)p(x)dx (9) 

J13.9 V 7T / 

where terms up to 0(a*) may be found in refs. [13], [I5|,[I5|. We use here the 
highest known order in a s because I 2 may be considered as an experimental 
input. The charm continuum contribution is then given by 



l cont 



h - h (10) 



Finally we have to take into account the two J/ty and charmonium 
resonances below the continuum threshold 

I res = 9tt(137.04) 2 {m x Y x p (m?) + m 2 T 2 p (m 2 2 )) (11) 

where the resonance masses and widths, are given by 

mi = 3.0969GeV 

r\ = (5.26 ± 0.37) x 10" 6 Ge^ (12) 
m 2 = 3.6860Ge1/ 
T 2 = (2.12 ± 0.18) x \Q- & GeV. 

The complete hadronic contribution to the sum rule is sum of these two 
integrals 

1 charm lcont ~\~ Ires (13) 
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5 Polynomial weight functions 

There is much freedom in the choice of the polynomial in the duality relation, 
Eq.(p. In this note we will use the polynomial to reduce the importance of 
the continuum contribution relative to the well established sub-threshold 
resonances. We impose the following conditions: 

The polynomial p(s) — 1 at s — mj^, it vanishes at the end of the 
integration range (s = 25Gev 2 ), and it is small in the continuum region. To 
be specific we choose an n-th order polynomial p n (s) = cto + ais + .. + a n s n 
and determine the coefficients by the constraints 

Pn(25) = 

PnW/*) = 1 ( 14 ) 

/ s k p k (s)ds = 0, k = 0, 1, .., n - 1 

J3.9 

For a 3-degree polynomial we obtain, for example 

p 3 (s) = 7.933 7875- 1.209157911s 

+6.015 360 076 x 1(TV (15) 
-9.792 537 729 x 10" V 

This result is plotted in Fig|| 




Figure 2: Polynomial fit vanishing at the continuum range (13.9 — 25)GeV 2 . 

It is obvious from the figure that the continuum data will almost cancel 
when integrated with this polynomial. 
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6 Results 

From the duality sum rule Eq. (|]) we obtain with the help of Eqs. (|8|-|T3|) and 
Eqs.flU) 



n 6 3 

' charm 



^J2J2 a k A ij( m c,^)( m f) J( k -hi) (16) 



fe=0 i=0 j=0 

which can be solved for m c . 

We present here predictions for simple duality, i.e. p(s) = 1, and for the 
3rd-degree polynomial above. The unknown m c in Eq. (|16|) is the running 
charm mass at a scale /!., which we fix to be /i = 5GeV, the relevant scale of 
the problem. 

For the coupling constant a s we take as an input its value at the mass of 
the tau lepton []I7J 

a s (m T ) = 0.345 ±0.020 (17) 



with m T = l.mGeV. After appropriate matching |L8| from 3 to 4 flavors 
this corresponds to 

a 4f (hGeV) = 0.224 ± 0.013 (18) 

For simple duality, we plot in Fig.^ the contribution of the QCD integral 
(lhs of Eq.(|l^)), at tree level and first and second order in the strong coupling, 
function of the charm quark mass m c . 




0.5 1 1.5 2 



Figure 3: QCD integral without polynomial weight as a function of mass. 
The three curves (from top to bottom at m = 2) represent tree, first and 
second order calculations in the strong coupling constant. 
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For the Data integral we have the continuum contribution I cont = lA.OQGeV 2 
and the resonance contribution I res = 12.80GeV 2 with their sum being the 
total charm data integral 



' charm 



26.86GeV 2 (19) 



Solving Eq.flTED for this value of charm data, the results for the mass 
of the charm quark at different orders in the perturbative expansion are: 
mf = OMQGeV, m« = 0.980GW and = 0.990GW. We see that the 
convergence of the QCD asymptotic expansion is extremely good. The main 
source of uncertainties come from the strong coupling constant and data. 
The result that we quote with this approach is 

m c (/i = 5GeV) = (0.99 ± 0M asymp ± 0.01 Qs ± 0.04 res ± 0.10 cont )GeV (20) 

The asymptotic uncertainty comes from the difference of two- and three-loop 
results used for the QCD correlators. The result of Eq. corresponds to an 
invariant mass |i~9f 

m c (m c ) = (1.40±0.11)GeV (21) 

In this final result the errors have been added quadratically. 

The contribution of the gluon condensate is completely negligible 
Our second approach consist in plug-in the 3rd degree polynomial in 
Eq.([il|) into Eq.([Tj]) in order to minimize the contribution from the charm 
continuum data, and therefore minimize the error involved in these data. As 
before we plot in Fig.|] the QCD integral as a function of m c for different 
orders in the perturbative expansion 

The contribution from the continuum data is now I cont = —0.17QGev 2 , 
whereas the resonance contribution, practically from J/ty resonance, is I res = 
9.341GeV 2 . The complete charm data contribution is then 



' charm 



9.165GeV 2 (22) 



With this value we solve again Eq. ( 16 ) at different orders in the perturbative 
expansion of a s with the results: = l.l53GeV , rrvp = 0.99lGeV , = 
0.927^6^. We see that, although we have eliminate the uncertainty coming 
from continuous data, the price we pay is that the asymptotic expansion does 
not converge so nicely as before, but still good enough to make a sensible 
prediction for the mass of the charm quark. We have 

m c (n = 5GeV) = (0.93 ± 0.06 a „ ± 0.015 as ± 0.0U res )GeV (23) 

This corresponds to an invariant mass |19| 

m c (m c ) = (1.34 ± 0.08)GeV (24) 
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m 



Figure 4: QCD integral with polynomial weight as a function of mass. The 
three curves (from top to bottom at m = 1) represent tree, first and second 
order calculations in the strong coupling constant. 

The influence of the gluon condensate is in this approach ~ 0.3% for {a s GG/7r) ~ 
.024GW 4 , still negligible. 

The result we find is perfectly compatible, within error-bars, with the one 
we found above by the simple duality approach. This agreement constitutes 
a non trivial confirmation of the duality ansatz. Averaging the results of 
both approaches, we finally find 

m c (m c ) = (1.37 ± 0.09)GeV. (25) 

Our value for m c (m c ) appears to be slightly bigger than the ones given by 
alternative QCD sum rule methods. With our normalization point the latter 
results read m c (m c ) = (1.29±O.O5)GeV0 and m c (m c ) = (1.28±0.06)GW@. 
Although these values agree, within error bars, with ours, it should be kept 
in mind that these authors use older and lower values of the QCD coupling 
constant, so that the error-bar quoted should really be larger. 

7 Conclusions 

In this letter we have analyzed the mass of the charm quark in the context of 
QCD finite energy sum rules. In the phenomenological side of the sum rule 
we use recent BESII e + e~ data, whereas in the theoretical side we employ 
the a large momentum expansion of QCD vector correlator function up to 
0(a 2 ) and 0{m 12 / q 12 ) . Two approaches are considered. The first one 
uses simple Cauchy sum rule for the correlator. The second one includes a 
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polynomial in the sum rule to minimize the contribution of the continuum 
data. The results from both approaches are nicely compatible with each 
other. Whereas with the first approach suffers from a substantial uncertainty 
arising from the continuum data, the second one shifts this uncertainty to 
QCD asymptotic expansion. The results allow a nice consistency check of 
QCD duality assumption. More precise results need either better data or 
further terms in QCD asymptotic expansion. 
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